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Abstract 

Dynamical systems theory provides powerful methods to extract 
effective macroscopic dynamics from complex systems with slow modes 
and fast modes. Here we derive and theoretically support a macro- 
scopic, spatially discrete, model for a class of stochastic reaction- 
diffusion partial differential equations with cubic nonlinearity. Di- 
viding space into overlapping finite elements, a special coupling con- 
dition between neighbouring elements preserves the self-adjoint dy- 
namics and controls interelement interactions. When the interelement 
coupling parameter is small, an averaging method and an asymptotic 
expansion of the slow modes show that the macroscopic discrete model 
will be a family of coupled stochastic ordinary differential equations 
which describe the evolution of the grid values. This modelling shows 
the importance of subgrid scale interaction between noise and spa- 
tial diffusion and provides a new rigourous approach to constructing 
semi-discrete approximations to stochastic reaction-diffusion partial 
differential equations. 

1 Introduction 



In modelling complex systems we often desire to model the effective macro- 
scopic dynamics. Some cases can be extracted from the full, or microscopic, 
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description by methods such as averaging, invariant manifold reduction and 
homogenization |H HS1 e.g.]. 

Stochastic partial differential equations (spdes) are widely studied in 
modeling, analyzing, simulating and predicting complex phenomena in many 
fields of nonlinear science [21 [3l [T5j |23j e.g.]. Recently, macroscopic reduction 
for dissipative spdes, with two widely separated timescales, has been studied 
by the dynamical systems theory of stochastic invariant manifolds [SJ EJ CEO 
e.g.] and by averaging methods [T9]. 1201 e.g.]. Moreover, invariant manifold 
theory also applies to generate a macroscopic discrete model of deterministic 
and stochastic dissipative pdes HOI El e.g.], the so-called holistic finite 
differences. Roberts [11] recently extended the approach to ensure macro- 
scopic discrete models preserve important self-adjoint properties of the fine 
scale dynamics for deterministic systems. 

Here we address the macroscopic discrete modelling of dissipative spdes 
and develop a novel rigorous approach. We consider reaction-diffusion in a 
one dimensional spatial domain driven by a noise which is white in time and 
with spatial structure. Let the non-dimensional spatial interval / = [0, L] 
with length L > 0, and let L 2 (I) be the Lebesgue space of square integrable 
real valued functions on I. Consider the following non-dimensional stochastic 
reaction-diffusion equation for a stochastic field u(x, t), of period L in space x, 

d t u = d xx u + a(u — u 3 ) + a d t W on J, (1) 
u(Q,t) = u(L,t) and u x (0, t) — u x (L, t), (2) 

where W(x,t) is an L 2 (I) valued Q- Wiener process defined on a complete 
probability space (ft, J 7 , P) which is detailed in the next section. 

The spatial domain I is divided into M elements and consequently a set of 
M fields are defined, one on each of these elements. In order to preserve the 
self-adjoint property of the linear operator defined on the elements we impose 
special interelement coupling conditions with a strength parametrised by 7. 
The system dynamics are expanded in the interelement coupling 7 so that 
based upon the case of weak coupling, that is, small 7 > , an asymptotic 
approximation and an averaging method derive a family of coupled stochastic 
differential equations (sdes) which describe the evolution of grid values; that 
is, the amplitude of the system on each element. These SDEs are a discrete 
stochastic model of the continuous space stochastic system ([I])-©. Our dis- 
crete SDE model highlights the macroscopic influence of 'subgrid' interactions 
between noise and spatial diffusion in the spde. 

The simplest conventional finite difference approximation of the SPDE (pQ) 
on a regular grid in x, say Xj = jh for some constant grid spacing h, is 

dUj = T2 (Efj+i - 2£/j + Uj-i) dt + a(Uj - Uf)dt + a dWj (3) 
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where Uj(t) is the grid value of the field u(x,t) at the grid points Xj, and 
similarly Wj(t) = W(Xj, t). However, our analysis herein recommends that a 
more accurate closure incorporates stochastic influences from the neighbour 
grids as in the Ito system of SDEs 

dUj = ^{U j+ i - 2Uj + Uj-x) dt + {aUj - aUf)dt + a dWj 

+ SyftUjdWj + ^ (dW j+1 - 2dWj + dWj-xj . (4) 

The second terms in the last line of the above SDE system reflect interaction 
between noise and spatial diffusion. Terms in Wj and a are due to the 
microscopic, subgrid scale, stochastic interactions discussed in Sections El 

In order to generate a macroscopic discrete model we divide the domain 
into finite overlapping elements and choose special coupling boundary con- 
ditions, Section El Such interelement coupling rules were first introduced by 
Roberts [TTJ to construct spatially discrete models of deterministic dynam- 
ics. One important property of this interelement coupling is the preservation 
of self-adjoint symmetry in the underlying spatial dynamics. Moreover, the 
strength of the coupling is parametrised by 7, < 7 < 1 : when the coupling 
parameter 7 is small, the coupling is weak and the system separates into 
'uninterestingly' decaying fast parts and the relevant slow parts, Section 
Then an averaging method [Tj5] derives a reduced model which describes the 
evolution of the overall amplitude of ([I])-© on the whole domain. Further, 
an analysis on eigenfunctions obtains a reduced model describing the evolu- 
tion of the local amplitude on each element. This model is the macroscopic 
discrete approximation to jl])-(j2]) expressed in (j4j). 

In this approach one difficulty is to construct from the original spatio- 
temporal noise W(x,t) a Wiener process Wj(x,t) on each element. A nat- 
ural method is expanding W(x,t) by the eigenfunctions of the linear op- 
erator £ 7 on each element, Section [21 which has analogues in the method 
of finite elements |16j . This construction shows that for small h and small 
coupling 7, the stochastic force on the first mode of £ 7 approximates the 
grid values W(Xj, t). Moreover, analysis of the eigensystem of £ 7 for small 7 
shows the slow parts of the system dominate: the fast parts converge to quasi- 
equlibrium with rate 1/h 2 . Then by this and the construction of IV 7 (x,t), 
the macroscopic discrete reduced model is proved to be consistent to the 
stochastic react ion- diffusion equations ([I])-© as the element size h — > . 

Another difficulty is that the linear operator £ 7 , defined in (fTU|) . varies 
with the coupling parameter 7. So the continuity of the linear oprator £ 7 
in 7 is needed. Section H] argues that the graph convergence of £ 7 as 7 — > 
ensures the continuity of eigenfunctions and eigenvalues in 7, then an asymp- 
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totic expansion of the first eigenfunction in 7 shows that the grid value is the 
amplitude of the system on the element. Then by averaging and the expan- 
sion of the first eigenfunction of £ 7 , the macroscopic reduced model is in the 
first eigenspace of £ 7 , which is varying with respect to 7. So last we project 
the reduced model to the first eigenspace of Co, the basic mode, to derive 
the macroscopic model. In this approach one interesting phenomenon is that 
the effect of noise in the subgrid scale fast modes is transmitted into the 
macroscopic slow modes by the projection. Numerical simulations confirm 
such transmittal [TO] . 

2 Overlapping finite elements and coupling 
boundary conditions 

This section divides the spatial domain I into M overlapping elements with 
grid spacing h. Let the jth element 

I j = [X j -h,X j + h] 

with grid points Xj — jh on each element Ij, j = 1, 2, . . . , M . Here for the 
periodic boundary condition we use the notation Xj± M = Xj . Let Uj(x) 
denote the field on the element Ij, j = 1, 2, . . . , M . Denote by f(u) = —u 3 . 
Then, modified from the spde ([Tj), consider the following system of spdes 
defined on elements Ij, j = 1, 2, . . . , M , 

d t u](x,t) = d xx u](x,t)+a7 2 u](x,t)+af(u](x,t))+ad t Wj(x,t) on Ij, (5) 

with the following interelement coupling conditions on the fields parametrised 
by 7, and with 7' + 7 = 1 , 

u](X J±1 ,t)= 1 'u](X 3 ,t)+ 1 u] ±1 (X J±1 ,t), j = 1,2,..., M, (6) 

and coupling of the first spatial derivative, denoted by subscript x, 

ul x (X~,t) - ulJXf, t) + Tul-ijXj, t) - Ju] +1 , x (Xj, t) 

-iuj^Xj-u t) + l'ul x {X j+1 , t) = , j = 1, 2, . . . , M , (7) 

with, to account for L-periodicity of solutions, 

u]± M (x ±L,t) = u](x, t) (8) 
u] ±M)X (x ±L,t) = u] x (x, t) . (9) 
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The coupling parameter 7 controls the flow of information between the two 
adjacent elements: when the coupling 7 = 0, adjacent elements are decou- 
pled; when 7 = 1, the system is full coupled and dSJ) (JTJ) is equivalent to the 
dynamics of the physical stochastic reaction-diffusion equation ([I])-©, see 
Section [3j 

The noise fields Wj(x, t) defined on each element Ij, j = 1, . . . , M , are 
infinite dimensional Wiener process which are detailed later from W(x,t). 

Related but different interelement coupling boundary conditions empow- 
ered an earlier exploration of the non-self-adjoint interaction between noise, 
nonlinear advection and spatial diffusion in discretely modelling the stochas- 
tic Burgers' equation [TO] . 

For our purposes, first we introduce a mathematical framework for sys- 
tem d5])-([7j). Let Hj = L 2 (Ij) be the set of all square integrable function 
on Ij and Vi = H x {Ij). Denote by U = U.f =1 Hj , V = Iif =l Vj and 

U a = Hf =l HJ , a > . 

Here H? denotes the usual Sobolev space W 2,a (Ij) [1 4 j . Define the inner 
product (•, ■) on % as the sum of the element integrals 



M 

(u,v) = 

3=1 



X- r x i+i 

Uj(x)vj(x) dx + / Uj(x)vj(x) dx 



for any u,v G H with u = (uj) and v = (vj) . And denote by || • || the 
product L 2 -norm on space "H. For any a G Z + denote the semi- norm on 'H a 

as 

M 

For the system ©-([Tj) we introduce the family of functional spaces 

H-y = {(Uj) G U : Uj(X j±1 , t) = iv.j{Xj, t) + Ju j±1 (X j±1 , t), 
j = 1, 2, ...,M, with L-periodicty (JB]) 



and the subspaces 

v 7 = {(a,),- vh- :»:.,.(.v ; ./)-^,. s .v/./) • 1U ) { jx r t) 

- 7uJ +M >(*;, t) - iu^iX^t) + Ay*;+i, t) = , 
j = 1, 2, . . . , M , with L-periodicty Q}. 

Then define the second order differential operator £ 7 : L>(£ 7 ) C V 7 — > H by 
C jU = for all u = { Uj ) G D{C y ). (10) 
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By a basic calculation [TT], — £ 7 is a self-adjoint second order operator. A 
direct calculation yields that for any u = (v,j) G D(C^) 

M 

(-£ 7 u,u) = ^ IKIli ^ ( n ) 

i=i 

which establishes the positivity of operator — £ 7 . Then there are coupling 
dependent eigenf unctions {(ej fc (x))}^ which form a standard orthonormal 
system in space % 7 and a sequence of real numbers < Ao(7) < ^1(7) < • • ■ 
such that 

= -A*(7)(cJ lJfc ), fc = 0,1,..., 0< 7 < 1 ■ (12) 

Moreover, £ 7 is an infinitesimal generator of a Cq semigroup; denote this 
semigroup by {<5> 7 (£)}<>o . By the positivity of — £ 7 we define (— £ 7 ) Q for any 
exponent a > by 



-£ 7 ) a w = [ x k(l)] a (u jjk e] tk 



k 



for u = (uj) = (J2k u j,k e ] t k) ■ Denote by "H 7 = D((— £ 7 ) a ) and define the 
semi-norm || • || a)7 in space by 

IMk 7 = |l( - £y) a/2u llo> «eV r 

By the same calculation as ffTTj) . for a G Z + U {0} , 

IM|a )7 = ||«||a, M G . (13) 

Given a complete probability space (f2, J 7 , {J-f}i>o, P), define the £ 2 (/) 
valued Q- Wiener process 

oo 

W r (z,*) = 5i) > /?fcM*)e*(aj) 

fe=0 

where {/3fc(£)}fc are mutually independent standard Brownian motions and 
{efc(x)}fc is a standard basis of L 2 (I) with eo(x) = yl/L and for > 1 

2 cog 2«p ^ k = 2m, 



2 gin 2»p ^ k = 2m — l. 



e k (x) = 

Moreover, assume that the Wiener process is sufficiently well-behaved that 

oo 

^%.<oo. (14) 

A:=0 
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Now we define the % 7 -valued Wiener process (Wj(x,t)) on all the over- 
lapping elements in the following series form 

oo 
1=0 

where {/3j,i}kLo are mutually independent standard Brownian motions on 
(Q, J 7 , P) and there are G R such that 

Here for each element j = 1, . . . , M 

yJ^MWiM = (^(x,t),ej,(o;))ej/x)/|| e jjo. 
Moreover, by the assumption (|14p 

oo 

J2MlH,i<B<oc, (16) 

1=0 

where the bound B is independent of the coupling parameter 7 G (0, 1]. 

Remark 1. The definition of W 1 (x, t) is similar to the definition of that in 
the finite element method [TBI I2H E5] ■ 

Now for fixed h > and for any T > , 

Wy(x,t) = (W7(x,t))eC 1 /\0,T;U*). 

and {W /7 (a;, t)}o< 7 <i is compact in space C(0, T;H°) for a < 2 for almost 
all w G . Then for almost all u G f2, the following limit 

M?(x,t) = (WAx,t)) = \im(Wj n (x,t)) 

is well defined in space C (0, T; %) for some 7„ — > 1 as n — > 00 . Moreover, 
by the coupling conditions we have almost surely 

Wj(X j±1 , t) = W j±1 (X j±h t), j = 1, 2, . . . , M , 

with notations W (0,t) = W M (X M ,t) and W M+1 (X M+1 ,t) = W^X^t) . 

Remark 2. By the analysis on eigenfunctions (ej fc (x))fc in Section m the 
above limit of W" y (x, t) in space C(0, T; %) is unique in the sense of distribu- 
tion for any sequence 7„ — > 1. Further, the distribution of Wj(x,t) coincides 
with that of Wj±i(x,t) on the common overlapping domain. 
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Finally, we explore the linear operator £ 7 as 7 — > , denoted by £ . 
Define 

£oM= (S%)' u =M* D fa)> (17) 
with the 'insulating' version of the coupling conditions (JO])-©, j = 1, . . . , M , 

u 3 {X ]±l )=u 3 {Xj)= Uj {Xf) ) 

u iA x t) ~ u 'jA x j) + u iA x i-i) - UjA x i+i) = • 

By a standard computation [II] . the spectrum of — £ > {^/cjfcLo > * s 

{0, 7r 2 //i 2 , 47r 2 //i 2 (triple), 9vr 2 //i 2 , 167r 2 //i 2 (triple), . . . , fcV/Zi 2 ,...}. 

Denote the corresponding orthonormal standard eigenmodes on each element 
by {e jtk }T=o > 3 = !> • • • > M > tnen e i,o(a;) = V and {e^- fc (x)} fc >i are 

f 1 ■ Tr(x-Xj) f 1 2it(x-X 1 ) 1 ■ 2tt(x-X,) 1 - 2tv\x~X 1 \\ 

{ 7K sm "A > 1 71 cos 71 sm > 7^ sm - fc ) > 

_ sm -4-^2, |^ cos -V 1 , ^ sm -V^> sm ~* J > 
...}. 

3 Limit system for full coupling 

Now we show that for full coupling, that is, as 7 — > 1 , equations ©-(JTj) gen- 
erates a model for the dynamics of the original physical stochastic reaction- 
diffusion equation ([I])-©. 

This is followed by a discussion similar to the case of Dirichlet boundary 
conditions [21] . here we just state the result and omit the detailed proof. 

Theorem 3. Assume bound f lT6|) and u 7 (0) e % 2 with ||w 7 (0)||2 < Co which 
is independent of '7 and u 7 (0) — > u° in % as 7 — > 1 . Tnen /or an?/ T > , u 1 
converges to u in distribution as 7 — )■ 1 m space C(0, T;%) where u = (v,/) 
solves 

d t Uj = d xx Uj + cra| + af(uj) + ad t Wj on Ij (18) 
TOi/i «(0) = u° and 

Uj(X j±h t) = u j±1 (X j±1 , t) , uj(x, t) = u j±M (x, t) . (19) 

To prove the above result needs some energy estimates on the solutions 
m 7 (x, t). By the definition of £ 7 , the SPDEs ©-([7]) takes the following ab- 
stract from 

du<(t) = [/> 7 (t) + a 7 V(t) + aF^it))} dt + adW^t) , 
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tt^(O) = (t»J(0)) , (20) 

where w 7 (t) = (uj{t)) and F{u') = (f(uj)). Writing the temporal depen- 
dence explicitly, in a mild sense we have 

ui{t) = S 7 (t)w 7 (0) + a [ S 7 (t - s) [7 2 it 7 (s) + F(w 7 (s))] ds 

Jo 

Jo 

Then by a standard semigroup approach |6j, for any w 7 (0) G H 7 and any 
T > , (JUD lias a unique mild solution u^{t) G C(0,T;ft 7 ) n L 2 (0,T; V 7 ) . 
Letting 

z i{f) = o [ <S 7 (t- s)dW' y (s) , 
Jo 

we have the following lemma. 

Lemma 4. Assume boundedness (TT61) . For any T > and q > i/iere zs a 
positive constant C q {T) such that 

E sup ||z 7 (t)||« 7 <C ? (T). 

0<t<T 

Proof. Since — £ 7 is positive and self-adjoint, 

||&y(t)|Lc<l, *>0. 
Then the result follows from the stochastic factorization formula [6]. □ 
Now define the difference w 1 = u 1 — z 1 , then 

dw^it) = [£ 7 u> 7 (t) + a 7 2 w 7 (t) + aFiu^dt , w 7 (0) = w 7 (0) . 

By the standard energy estimate to stochastic reaction-diffusion equations 
with more general nonlinearity [19J we have 

Lemma 5. Assume m 7 (0) G "H 7 , i/ien /or any T > 0, there is a positive 
constant Ct > such that for any p G Z + 

E||M 7 ||c(0,T;« 7 )niP(0,T;V 7 ) + E|| <%« 7 || L2( 0) T;«-i) < Ct(1 + ll« 7 (0)||o) . 

We show that {"D(-u 7 )} 7 , the distribution of u 1 in space C(0, T;H), is 
tight. For this we need the following lemma by Simon [T3] . 
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Lemma 6. Assume E, E and E\ are Banach spaces such that E\ d E , 
the interpolation space (E , Ei) e ^ C E with 9 G (0,1) and E C E with 
C and (s denoting continuous and compact embedding respectively. Suppose 
Po,Pi G [1, oo] and T > , such that 

X is a bounded set in L P1 (0,T; E\) 

and 

dX := {dv.ve X} is a bounded set in L Po (0, T; E ). 
Here d denotes the distributional derivative. If 1 — 9 > 1/pe with 

1 1 — 6> 9 

Pe Po Pi ' 

then X is relatively compact in C(0,T; E). 

By the above lemma, and noticing the relation (I13|). we have the following 
theorem. 

Theorem 7. Assume ffT6]) and u 7 (0) G H 7 wift ||u 7 (0)||o < Co which is 
independent 0/7. For any T > , T>(u 7 ) zs tzo/ii m C(0, T;7i) . 

Similarly if w 7 (0) G "H^ with ||w 7 (0) || 2)7 < C which is independent of 7, 
for any T > there is a positive constant Ct > such that 

E sup ||M 7 (t)|| 2 , 7 < C T - 

0<t<T 

Then by the embedding of H 2 (I) C C^i") 
E 



du] , 

<C T , 0<t<T. (21) 



T^fPW) 

By the above estimates we can treat the boundary value in passing to the 
limit 7 — > 1 of full coupling. 

By Theorem [TJ for any k > there is a compact set K K C C(0,T;H) 
such that 

p{m 7 g iq > 1 - « . 

Then there is a function w G C(0, T; %) and a subsequence 7„ — )■ 1 as n — > 00 , 
such that in probability 

■u 7 " — > u as n — > 00 . 
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Now we determine the equation solved by the limit u. Define a test function 
ip E C °°(0,L) and define 

tfj = (pli, ■ 

Then by the boundary conditions we have in the variational form for the 
system ©-© 



(u* (f) >V >) 



(u*(0), <p) + / (C Jn u^(s) + atfu*(s), <p) ds 
Jo 

+ / (aF(u' yn (s)),ip)ds+ [ (adW ln {s),(p) 
Jo Jo 



■ JO JXj-! 



dx dx 



dx ds 



t rXj +1 

Jxj-i 
t r X 



dWj n (s)ipjdx 



/ / Uj n (s)ipj dx ds + a \^ / / f{uJ l {s))Lpjdxds 

a Jo Jxj-x ' , Jo Jxa_x 
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n 



o L 



du^X^s) duf +1 (X v s) 



dx 



dx 



duf-iXj.ua) du]-(X j+u a) 



dx 



dx 



¥j{Xj)ds 
(fj(Xj) ds . 



(22) 



Then by estimate fl2Tl) . letting n — > oo , that is 7^ — > , the last two terms 
disappear. Notice that f(u] n ) — » /(it?) weakly in space L 2 (0,T;L 2 ) and by 
the assumption on W 7 ™ we have, by passing to the limit n — > 00 , 



' ( %(s) d<pj 

dx dx 



dx ds 



t r X 



3+1 



Jx 



(uj(s) — u?(s))(p dx ds 



+ / (<TdW{s),<p) 



(23) 



with = (Wj(t)) which is well defined by Remark [2] Then a density 

argument yields that u = (v,j) solves the following stochastic equations 



d t Uj = [d xx Uj + aiij — a(uj) 3 ] dt + a d t Wj on Ij, 
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with coupling boundary conditions 

Uj(Xj±i, t) = Uj±i{X j±1 , t) , Uj(x, t) = u j±M (x ±L,t). 

Remark 8. By the boundary condition (|T9l) and Remark[2], the distributions 
of Uj and Uj+\ in space (7(0, T; L?(Xj,Xj + h)) coincides. 

Now define 

u(x, t) = Uj(x, t), x G [Xj, Xj + h] 
and an L 2 (/)-valued Wiener process W(x,t) as 

W(x,t) = Wj (x, t), xe [Xj , Xj + h] . 

Then u(x,t) solves the stochastic reaction-diffusion equation (DQ) — (pD with 
the noise term W(t) replaced by W(t) without changing the distribution. 
So recovers the original system ([I])-©, in distribution, in the limit 

of full coupling, as 7 —7- 1 . 

4 Amplitudes on the elements 

We derive a discrete macroscopic approximation to the system of SPDEs (^j- 
(J7J) based upon small coupling parameter 7 > . By the analysis on opera- 
tor £ in section^ for 7 = the dominant mode is (e j0 ), so by hyperbolicity 
we expect that for small 7 > 0, the dominant mode is (ej ). This is followed 
by the analysis on the continuity of {(ej fc )}fc and Xk{j) on coupling param- 
eter 7. Further the asymptotic expansion for (ej ) in 7 shows that the grid 
value approximates the amplitude on each element. 

For this we study the continuity properties of £ 7 as coupling 7 — > . We 
use variational convergence for operators pQ. For any subsequence 7 n with 
7„ — > as n — > 00 , we introduce the G-convergence for £ 7n . 

Definition 9 (G-convergence). Operator C ln is said to be graph- convergent 
(G- convergent) to Cq as n — >■ 00 if for every (u,v) with v = CqU, there exists 
a sequence (u n ,v n ) with v n = C ln u n such that u n — >■ u strongly in V and 
v 11 — > v strongly in V* , the dual space ofV. 

Now for any u = (uj(x)) G Vo, denote by v — CqU . First we choose 
bounded set {v n } C "H 7n such that v n — > v in the dual space V*. Then solve 
the following equation 

C Jn u n = v n . (24) 

By the relation fTl3|) . {u n } is bounded in 1-L 2 which yields that {u n } is compact 
in V. Then there is a subsequence, which we still denote by {u n }, that 
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converges to u in V as n — )■ oo . Multiplying testing function ip G C^°(0, L) 
on both sides of (T23J) and passing to the limit n — > oo , we have 

Cqu = v (25) 

which yields that u = u by the uniqueness of the solution to ( I25p . Then we 
have £ 7n is G-convergent to £ • 

Now we draw the following result on the continuity of eigenvalues and 
eigenfunctions in coupling 7 [2T] . 

Theorem 10. 

lim A fc ( 7 ) = A fc = -k 2 7r 2 /h 2 , k = 0, 1, 2, . . . . 

Lei mfc fre the multiplicity of \k(l), the sequence of sub spaces of dimen- 
sion mk generated by ( (ej£ ),..., (ej^ mfc )) converges in % to the eigenspace 
of Cq corresponding to A& = —k 2r n 2 /h 2 . 

Remark 11. Such convergence of the eigenspaces is called Mosco- conver- 
gence pp. In the sense of this convergence, (ej^) may converges to (e^' fc ) 
instead of (e l j k ) , 1 < I ^ V < . 

By the above result for small coupling 7 > , we write for each element 
j = l,...,M, 

e lo( x ) = e Jfi(x) + ejo(x) 

with 

^Jo( x ) — as 7 . 

Next we show that ej (a;) has an expansion to 7 2 terms in the coupling 
parameters: 

d 2 

— q {x) = x {j)q {x) + A ( 7 )/v / 2l, 

with boundary condition (|S])-([7J) • Now if we account for the coupling and 
boundary condition by 

eJ,o(z) =eJ,oW+?,o(^) 



+ 



2 [eJ i0 (X,) - eJ_ 1)0 (X i _ 1 )] (* - X,) , X,-_ a < x < X j , 
. g [e]+ifi{Xj+i) ~ *U X ij\ ( x ~ X ^ > Xj<x< X, +1 . 



Then e^x) solves 

9 2 



(x) = Ao( 7 )e7, (x) + A ( 7 )/V2/i (26) 



dx 2 3 
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with boundary condition 

Ha(X j±1 ) = e jfi (Xj) = 0. 

But, by (j7j) we also have coupling in the derivatives: 

d x e] >0 (Xr) - d x e] fi (X+) + -<),,] L0 (.Y,-) - 
'U-,,(.V ; :) • -U^. (.V,.,) 

= 2^ [e] +lj0 (X j+1 ) - 2el (Xj) + e]_ lfl {X 3 ^)] , j = 1, 2, . . . , M , 

Then by the above boundary condition, equation fl26|) and the fact eJo( a; ) - ► 
as 7 — )■ , we have 

A o (7) = 0(7 2 ) as 7^0 

which implies that ej (x) have an expansion to 7 2 terms in the coupling 
parameter. 

Assume we have the following asymptotic expansion for each element, 
j = l,...,M, 

e], (x) = el (Xj) + 7 i^(x) + j'F^x) + Fj fl (x) (27) 

where -^3(2) = C( 7 3 ) • By Ao(0) = and the coupling boundary condi- 
tion ([6])-([7j) , Fj k , k = 1, 2 , are fcth order polynomial in x . Then also by the 
boundary condition -([7]) we have 

\ ™ J 7 " o( ^ (x-X,), X^x^X^, 

and 

j ' \A j (x-X j )(x-X j+1 ), X 3 <x<X J+1 , 

, ^oP^i) - ^IoW + e] +lfi (X j+1 ) 
3 ~ 2h* 
The above asymptotic expansion shows that for small coupling 7 > , the 
first mode is dominating and the grid value u^(Xj,t) approximates the am- 
plitude of the field w 7 (x,t) on the element Ij. 



with 
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5 Macroscopic models for small coupling 

By the asymptotic expansion in the previous section H] we derive a discrete 
macroscopic approximation model to (J5]) (JTJ) for small coupling 7 > . For 
this we first apply an averaging method to reduce (J5j)— (JT]) onto the slow 
mode (ej ). 

We split (it 7 ) into slow part and fast part. Define map Pj on H 7 to H. 



n 7 K 7 ) = (K,eJ i0 )eJ j0 (a;)/||eJ i0 ||g) and P 7 



T - P 7 
1 M) ) 



where / is the identity operator on "H 7 . And for no coupling, 7 = 0, write 
P = Pq and Pi = P® . Then denote by (uj(x, t)) the solution to system 
and make the following expansion 



k=0 



Now define the slow part and fast part, respectively, 

U\x,t) = (al(t)e] }0 (x)) and V\x,t) = (u](x,t)) - U\x,t). 
Then we have that these satisfy the coupled spdes 



dU 1 
dV~< 



C~ f W + «7 2 f/ 7 + aP^F{U\ V" 1 ) dt + a^dB, 



'0 ) 



£ 7 V 7 + a 7 V 7 + aP?F(U\ V 1 ) dt + a^dB} 



where 



00 

Bl{t) = (y^^oW) and Bj(t)=(j2^k^At)el 



ix, 



fc=i 



(30) 
(31) 



By the analysis of Section HJ ^0(7) = C(7 2 ) as 7 — > 0, then for small 
coupling 7 > , (j3~0~I) - (l3"Tl) have completely separated time scales. Thus an 
averaging approach applies to derive a macroscopic reduced system over the 
time scale 7~ 2 T for any T > [T5]. For this introduce a slow time scale 
t' = 7 2 t and small fields 



XP{t) = iW^H) and V^{t) = 7V A7 ( 7 2 t), 



(32) 



then on the slow time scale t' 



dW{H^ 



7" 2 £ 7 f/ 7 + aW + aPjF(W, V~<) dt' + a^dB^t') 
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dV\t') = 1 - 2 £ 7 V' + aV' Y + aP?F(U\V^) dt' + a^dB^t'). 

Here B^(t') = 7-Bo(7~ 2 <0 and Bj(t') = jBj(j~H') are Wiener processes 
with the same distributions as those of B^it') and Bj(t'), respectively, due 
to the scaling properties of the Wiener process. 

Let fp{t') = (^Jo(t')) £ -Pi 7 "H 7 be the unique stationary solution of the 
following linear equation 

dfp{t') = 7 - 2 £ 7 rf (t') dt' + (X7- 1 dBj(t'). (33) 

Then by an energy estimate and almost the same discussion to that by Wang 
and Roberts [19] , for any T > , there is a positive constant Ct such that 



sup E||Vr(0 - 7f(t')||o < 1 2 C T (E||^(0)||o + 11(^(0))© . (34) 

0<t'<T 

We have no explicit expressions of (ej fc (x)) , k > 1 , so for our purpose we 
derive another approximation for V r7 (t / ). 

Lemma 12. Assume bound ({TBI . For any T > 0, there is positive con- 
stant Ct such that 

sup E\\V^(t') - f}(t')\\ < 7 C T (E||^(0)||o+E||^(0)||o + ||(tt,- (0))||S) , 

0<t'<T 

where fj(t') = (fjjo(t')) is the unique stationary solution of the following linear 
equation 

dfj(t') = ^ 2 C fj(t') dt' + dB^t') 

and 

00 

Bi{t') = Y,{l^Ml- 2 t>M) 

k=l 

with distributions independent of coupling parameter 7. 

Proof. Expand fp{t) and fj(t) by the eigenfunctions of £ 7 and Co respectively 
as 

00 00 
ff(x, t') = Vl(x, t') and f)(x, t') = ^ Vk(x, t'), 

k=l k=l 

where 

rj k (x,t') = (fj(t'),(e j>k (x)))(e jjk (x))/M. 
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Then for k > 1 



and 



with 



dfg{t') = -^XMvUt')dt' + a -dBl k (t') 



dfjk(t') = -^X k fjk(t')dt' + -dB 1>h (t') 



BU*) = 7(^/3,, fc (7-V)eJ ifc (x)), 

Using the Ito formula and the stationary property of fjj and fjk , there is a 
positive constant C such that for k > 1 

E\\fjKt')\\o<C and E\\fj k (t')\\ <C. (35) 

Define the difference z^it') = fp k {t') — fjk(t') which solves 

dzl{t') = -1 [A*^(0 + (A*( 7 ) - hW k (t')\ dt' + °-d [^(0 - S^f) 
and hence 

r-t' 



7 Jo 

7 Jo 



Then by the analysis of Section @] on X^j), k >1 , assumption (|T6"|) and the 
estimates (1331) we have 

sup E||7f(0 -rKOUo <lC T E (||^(0) ||o+ ll^(O)Ho). 

0<t'<T 

Thus by (f3"l|) . the proof is complete. □ 
Now by (13"iZl) and the above result, for V 1 on the original time scale, 

sup E||^(t)- 7 ^(7 2 t)Ho < 1 2 C T (E||^(0)||o + E||^(0)|| + ||(%)(0))!lo) • 

0<t<7" 2 T 

(36) 
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Moreover, 

oo 1 

Ef(t) = a 2 ^— (g^ )fc (x)). (37) 
k=l k 

For IP we follow an averaging approach [19] which yields the following 
averaged equation 

dtP (If) = l^C^Wit') + aW(if) + aP^F(U^(t'))] dt' + ap 1 dW Q (t') (38) 
where U' 1 = (wj ) and 

F(-)=E[F(-,r/H)]. 
By the definition of / and that fj is Gaussian with zero mean, we have 



Moreover, by a deviation argument [333 |20], stochastic effects in these 
subgrid scale fast modes are fed into the slow modes by the nonlinear inter- 
action. So we have the following averaged equation plus deviation 

dlP (If) = [p 2 C^W + aW + aP^F(P 1 )]dt' 

+ ap 1 dBj(i') + cryy/Q(W)dP(t') (39) 

with, for fixed U 1 ', 

POO 

Q(W) = 2E / P 7 [F(W + fj(s)) - F(W)] 
Jo 

® P 7 [F(W + jj(0)) - F(W)] ds , 

and /3(t) = (j3j(t)) is an M dimensional standard Brownian motion. For 
any T > and any k > , there is a positive constant C K) t such that 

P {sup < t ,< r |£T(t') - t/ 7 (t')l < l 1+ C k , T } > 1 - k . 

Then for the original system fl5]), by (|32|) . we have the following reduced 
equation 

dW (t) = [P,U^(t) + a-fU^t) + aP^F(U^(t))} dt 

+ *jdBZ(t) + a 7 3 V 'Q{P l U2(t)) d$(t) , (40) 

where j3(t) = ($j(t)) = •y~ 1 (3('j 2 t) is the scaled M dimensional standard 
Browian motion. Then for any k > , there is a constant C Kj t > such that 

p { S up Q < t < 7 - 2T |cr(t) - uj(t)\ < 7 2+ a,T} > i - « . 
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Now having the reduced system (|40p . we use the approximation to the 
amplitude on each element to derive a further approximate model. Write 
U2{x,t) = (uj (x,t)) and define 

(Uj(t)) = al{t) (eT (Z,)) . 
Then by (l2"Tj) for (Uj (x,t)) we have the following asymptotic expansion 

«£„(*, t) = a2(t)eJ i0 (X,)+ 7 a2(t)^ 1 (x)+ 7 2 a2(t)i^ 2 (x) + O( 7 3 ) 

(UjW + MUjW-U^ityix-Xj) 

+ i(Uj-i(t) - 2Ui(t) + U j+1 (t))(x - Xj)(x - Xj^) 
+ G( 7 3 ), X,^<x<X 3 , 
U 3 {t) + l{U ]+l {t)-U 3 {t)){x-X 3 ) 
+ 2&(^-iW - 2U M + U j+ i(t))(x - Xj){x - x j+1 
+ 0( 7 3 ), Xj<x<X J+1 , 

Putting OH]) into (HQJ) yields 

, 2 



(41) 



dUjit) 



+ 0( 7 3 ,a 2 ), 



(42) 



where 



a, 



E [(^(s) - E?g , e,, )(^ (0) - Er/ 2 , e i)0 }] ds . 



Here we use the approximation of Gj<a{x) to ej (x) for small 7 . 

Notice that system (j4"2"j) is not a complete discrete model because the 
noise terms are still described on the mode (ej (x)). In order to give a 
discrete approximating model for small coupling 7 , we explore the evolution 
of the amplitude of the basic mode (ejfl{x))- For this we project (Uj) onto 
the basic space E spanned by (e J - i o(x)). However, the fast modes V 1 have a 
nonzero projection in basic space E ; because of the complicated expression 
for (ej k (x)), we choose to project 777, which approximates V 1 up to error 
of 0( 7 2 ). For this we first project ■yfj onto (ej (x)), then project to E . 

Notice that for small coupling 7, -fj(t') behaves as a noise process. By a 
martingale approach [5] [T9], [22] we have the following lemma. 
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Lemma 13. Assume bound (|T6|) . Then - J * fj(s) rfs converges in distribution 



fe=i v fc 

where (flj,k(t')) , k = 1,2 ... , are mutually independent standard M dimen- 
sional Brownian motion in time scale t' . 



Then on the right-hand sides of (1421) there are additional noise forcing 
terms when projected to the basic space Eq: namely 



oo h 

El 
V h 

k=l v fc 



jA x )> e lo( x )) d Pj,k{t) ej >0 (x) 



where 0j,k(t)) — (l~' l ^j,k(l 2 't)) , k = 1,2,..., are mutually independent 
scalar standard Brownian motions. By the expansion of (ej (x)) , (|27[) and 
the expression of (ej t k(x)) , 



(e Jifc (x),eJ i0 (x)) = 1 (e Jtk (x),Fl 1 (x)} + 0( 1 2 ). 



Then define 



fe=i 



and for j = 1 , 2 , . . . , M , i = -1 , , 1 

= VQj/3j(t)ej-ifi(Xj_i). 

Projecting the above system onto the basic space Eq and by (I37p we then have 
the following macroscopic discrete approximation model to the spde ©-(JT]) 
for small coupling 7 > 

dUj(t) = j- {Uj^it) - 2U j {t) + U j+1 {t)) dt + ^ajUjit) - aUf(t) 
+ <T-ydB jt0 (t) + a^dBj^t) + 3j 2 y/2Uj(t)dB jy (t) 



(T7 



(dB jtl (t) - 2dB j)0 (t) + dB jr 



+ ^{dBl 1 (t)-2dBl (t) + dBl_ 1 
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+ -^-l 3 Uj(t) (dB jtl (t) - 2dB j)0 (t) + dB jt _^j 

+ 0( 7 3 ,a 2 ). (43) 

Furthermore by (1251) . Fj-y = 0{^) . Then a truncation of (|43p to errors 
0( 7 3 , a 2 ) and evaluating at full coupling 7 = 1 yields the following macro- 
scopic discrete system of SDEs, j = 1, . . . , M , 

dUj(t) « ^{U j _ 1 (t)-2U j (t) + U j+1 (t))dt+(a j U j (t)-aUf(t))dt 
+ adB j)0 (t) + 3V2Uj(t)dB jfi (t) 

+ °- (dB hl (t) - 2dB ji0 (t) + dB^ . (44) 

This system of SDEs reduces to the system (T4]) discussed in the Introduction 
as a discrete model of the reaction-diffusion spde (Tj[|). 



6 Consistency of macroscopic discrete model 



Next we study the consistency of the macroscopic discrete model f l44j) by the 
definition of the Wiener processes W 1 {x, t) in Section [21 as h, the size of each 
element, converges to zero. 

By the properties of ej^(x) and ejfc(x) for small coupling 7 > in section |2j 
and section H] respectively, we have 



leTollo 



00 

= Yl VVkPkit) (e h (x), e] >0 {x)) e} i0 (x)/||ej j0 || 

fc=0 

00 

= £ yfaPk® (e k (x), el Q (X 3 ) + jF jtl (x) + i 2 F h2 {x) + ift(z)) ej )0 (x)/||ej >0 || 

fc=0 

00 

= <MPk(t) (e k (x), e jfi (Xj) + >yF jA (x) + i 2 F h2 {x) + F? 3 (x)) e jfi (x) + 0(j) 
00 

OO 

+ E 



k=0 



'2g 2 m-i - ...LI 



m=l 



L 1 

2mir 2h 



2mix 2miT 
sin — - — Xj+i — sm — - — 



m=l 

+ o(h n ) 



2mix 2h 



2mix 2mn 
cos — - — Xj-i — cos — - — X. 
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^vW*HM + 0(M- ( 45 ) 



k=0 

And similarly we have for I = 1,2, 



oo 

y/QkPk(t) (e k (x), e] ;l (x)) ej/x)/||ejj 



k=l 
oo 



Y V^kPkit) (e k (x), e j: i(x)} e j: i(x) + 0(7) 



k=l 



= 0(h,*y). (46) 

Theorem 14. The macroscopic discrete model (|44p zs consistent to the orig- 
inal stochastic reaction- diffusion equation (CO)-©. 

Proof. Following the exact same discussion by Roberts et al. [12] for deter- 
ministic systems, we just consider the stochastic terms and &j . 

By the definition of &j, since = 0[hj and noticing that = k 2 n 2 /h 2 , 
then we have 

oij = a + 0(h 2 ), h^O. 

For the stochastic terms, notice that fjk converges to zero in mean square 
with speed exp{— X k } = exp{—k 2 ii 2 /h 2 } as h — Y . Then since q^ x = 0(hj , 
we have Qj = 0(h 2 ^j as h — >■ . Then by the definition of W /7 (x, t) and analy- 
sis (SSD — (SHD the macroscopic discrete model (1441) is consistent up error O (h) 
to the following stochastic reaction-diffusion equation 

d t u(x,t) = d xx u(x,t) + a(u(x,t) — u 3 (x,t)) + o~B(x,t), x <G [0, L] 
u(0,t) = u(L,t) t>0 

where B(x, t) distributes the same as W(x, t) . This completes the proof. □ 



7 Conclusion 

Stochastic averaging is an effective method to extract macroscopic dynamics 
from SPDEs with separated time scale [T9J [2Q] - Here by applying the stochas- 
tic averaging and dividing the spatial domain into overlapping finite sized 
elements with special interelement coupling boundary conditions flS])-([7|), we 
derive a macroscopic discrete model f|44|) for stochastic reaction-diffusion par- 
tial differential equations (PQ) with periodic boundary conditions. The most 
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important property of such interelement coupling boundary conditions is 
preserving the self- adjoint symmetry which is often so important in applica- 
tion [12]. Furthermore, by the choice of stochastic forcing on each element, 
this coupling boundary conditions also assures the consistency for vanishing 
element size, section [6j 

Moreover, the final discrete model (1441) . which is different from the usual 
finite difference approximation model ([3]), shows the importance of the sub- 
grid scale interaction between noise and spatial diffusion and provides a new 
rigorous approach to constructing semi-discrete approximations to stochastic 
react ion- diffusion spdes. 
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